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y—i ' ABSTRACT 

o : 

( Context. Observations in the cosmological domain are heavily dependent on the validity of the cosmic distance-duality 

(DD) relation, Dl(z)(1 + z) 2 / Da(z) = 1, an exact result required by the Etherington reciprocity theorem where Dl(z) 

■ and Da(z) are, respectively, the luminosity and angular diameter distances. In the limit of very small redshifts Da(z) = 
QJ • Dl(z) and this ratio is trivially satisfied. Measurements of Sunyaev-Zeldovich effect (SZE) and X-rays combined with 

the DD relation have been used to determine Da(z) from galaxy clusters. This combination offers the possibility of 
, testing the validity of the DD relation, as well as determining which physical processes occur in galaxy clusters via their 

■ shapes. 

Aims. We use WMAP (7 years) results by fixing the conventional ACDM model to verify the consistence between the 
validity of DD relation and different assumptions about galaxy cluster geometries usually adopted in the literature. 
Methods. We assume that r\ is a function of the redshift parametrized by two different relations: rj(z) = 1 + rjoz, and 
r/(z) = 1 + r)oz/{l + z), where r/o is a constant parameter quantifying the possible departure from the strict validity of 
the DD relation. In order to determine the probability density function (PDF) of rjo, we consider the angular diameter 
distances from galaxy clusters recently studied by two different groups by assuming elliptical (isothermal) and spherical 
(non-isothermal) /3 models. The strict validity of the DD relation will occur only if the maximum value of rjo PDF is 
I ■ centered on r)o = 0. 

Results. It was found that the elliptical j3 model is in good agreement with the data, showing no violation of the DD 
' relation (PDF peaked close to 770 = at la), while the spherical (non-isothermal) one is only marginally compatible at 

3(7. 

Conclusions. The present results derived by combining the SZE and X-ray surface brightness data from galaxy clusters 
with the latest WMAP results (7-years) favors the elliptical geometry for galaxy clusters. It is remarkable that a local 
_ 1 ' property like the geometry of galaxy clusters might be constrained by a global argument provided by the cosmic DD 

relation. 
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1. Introduction 

■ The most useful distances in cosmology are the luminosity distance, Dl(z), and the angular-diameter distance, Da(z). 
' The expressions of both distances depend on the world models, but the relationship between them, namely 

X 



is valid for arbitrary spacetimes, a result usually referred to as distance-duality (DD) relation. 

The above expression can easily be deduced in the context of Friedmann-Robertson- Walker (FRW) cosmologies 
(Weinberg 1972). However, as originally proven by Etherington (1933), it depends neither on Einstein field equations nor 
the nature of matter content filling the spacetime. The proof depends crucially on photon conservation (transparency of 
the cosmic medium) and that sources and observers are linked by null geodesies in a Riemannian spacetime. 

The DD relation plays an essential role ranging from gravitational lensing studies to analyses of the cosmic microwave 
blackbody radiation (CMBR) observations, as well as for galaxy and galaxy cluster observations (Schneidder, Ehlers & 
Falco 1999; Komatsu et al. 2011; Lima, Cunha & Alcaniz 2003; Cunha, Marassi & Lima 2007; Ribeiro 1992, 2005; Ribeiro 
& Stoeger 2003). Indeed, any observational deviation from Eq. (Q]) would be a theoretical catastrophe thereby igniting a 
major crises in observational cosmology (Ellis 1971, 2007). 

Although taken for granted in virtually all analyses in cosmology, the DD relation is in principle testable by means 
of astronomical observations. One may assume a redshift dependence of the form 

^(l + z)- 2 =?? (z), (2) 
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where r)(z) quantifies a possible epoch-dependent departure from the standard photon conserving scenario (j] = 1). 

Basset & Kuns (2004) used both supernovae la data as measurements of the luminosity distance Dl and the estimated 
Da of FRIIb radio galaxies (Daly & Djorgovski 2003) and ultra compact radio sources (Gurvitz 1994, 1999; Lima & 
Alcaniz 2000, 2002; Santos & Lima 2008) in adopting this kind of approach to test possible new physics. Any source 
of attenuation ("gray"' intergalactic dust) or exotic photon interaction must violate the DD relation (More et al. 2009, 
Avgoustidis et al. 2010), thereby providing new consistency checks of cosmological models. 

On the other hand, observations of the Sunyaev-Zeldovich effect (SZE) from galaxy clusters are becoming a powerful 
tool in cosmology (Sunyaev & Zeldovich 1972; Cavaliere & Fusco-Fermiano 1978, De Filippis et al. 2005, Cunha et al. 
2007; Nord et al. 2009; Basu et al. 2010). The combination of SZE and X-ray provides the Da(z) of galaxy clusters, 
hence can be used to constrain some cosmological parameters. However, Uzan, Aghanim & Mellier (2004) argued that 
this technique is strongly dependent on the validity of the DD relation. When the DD relation does not hold (77 ^ 1), 
the observationally determined angular distance must be replaced by the value (in the Uzan et al. (2004) notation the 
correcting term is rj~ 2 ) 

Di ata (z) = D A (z) V \ (3) 

this quantity reduces to the conventional angular diameter distance (assuming cosmic transparency) only when the DD 
relation is strictly valid (7/ = 1). To quantify the 77 parameter, Uzan et al. (2004) obtained Da{z) as given by the cosmic 
concordance model (Spergel et al. 2003), whereas for D A iata (z) they considered 18 angular diameters from the Reese et 
al. (2002) galaxy cluster sample for which a spherically symmetric cluster geometry was assumed. By assuming 77 to 
be constant, their statistical analysis provided rj = 0.911q'q4 (la) and is therefore only marginally consistent with the 
standard result, 77 = 1. 

De Bernardis, Giusarma & Melchiorri (2006) also searched for deviations from the DD relation by using the D^ ata (z) 
from galaxy clusters provided by the sample of Bonamente et al. (2006). They found a non violation of the DD relation 
in the framework of the cosmic concordance ACDM model. Avgoustidis et al. (2010) adopted a DD relation of the form 
Dl = Da(1 + z) 2+e to constrain the cosmic opacity by combining the SN Type la data compilation of Kowalski et al. 
(2008) with the latest measurements of the Hubble expansion at redshifts in the range < z < 2 (Stern et al. 2009). By 
working in the context of a fiat ACDM model, they found e = -0.04+[j;^ (2a). 

In the past few years, many studies based on Chandra and XMM observations have shown that in general galaxy 
clusters exhibit elliptical surface brightness maps. Simulations have also predicted that dark matter halos show axis 
ratios typically of the order of « 0.8 (Wang & White 2009), thereby disproving the spherical geometry assumption 
usually adopted (Reiprich & Boringer 2002; Bonamente et al. 2006, Shang, Haiman & Verdi 2009). In this line, the first 
determination of the intrinsic three-dimensional (3D) shapes of galaxy clusters was presented by Morandi, Pedersen & 
Limousin (2010) by combining X-ray, weak-lensing, and strong-lensing observations. They studied the galaxy cluster 
MACS J1423. 8+2404 and found a tri-axial galaxy cluster geometry with DM halo axial ratios 1.53 ±0.15 and 1.44 ±0.07 
on the plane of the sky and along the line of sight, respectively. 

In this letter, we take the validity of the DD relation for granted to access the galaxy cluster morphology. The values 
of Da(z) are obtained from the WMAP (7 years) results by fixing the conventional flat ACDM model whereas the 
observational measurements D^ ata (z) are the angular diameter distances from galaxy clusters obtained via SZE plus 
X-ray techniques. These samples differ in terms of the assumptions concerning the possible cluster geometries of elliptical 
and spherical models. Our analysis is based on two parametric representations of r/(z) defined by Eq.[2](or Eq.[3]), namely 

1. J](z) = 1 + rjoz and II. 77(2) = 1 + r\Qzj(\ + z). 

The first expression is a continuous and smooth one-parameter linear expansion, whereas the second one includes a 
possible epoch- dependent correction, which avoids the divergence at extremely high z. These deformations of the DD 
relation effectively parametrize our ignorance of the underlying process responsible for its possible violation. However, 
we emphasize that these expressions are very simple and have several advantages such as a manageable one-dimensional 
phase space and a good sensitivity to observational data. Clearly, the second parametrization can also be rewritten as 
r\(z) = 1 + 7/o(l — a), where a(z) = (1 + z)^ 1 is the cosmic scale factor. This represents an improvement with respect 
to the linear parametrization, since the DD relation becomes bounded regardless of the redshift values. It will become 
more useful once higher redshift clusters data become available. 

The above parametrizations are clearly inspired by similar expressions for the o;(z)-equation of state parameter of 
dark energy models (Padmanabhan & Choudury 2003; Linder 2003; Cunha, Marassi & Santos 2007; Silva, Alcaniz & 
Lima 2007). In the limit of extremely low redshifts (z « 1), we have 77 = 1 and Dl = Da as should be expected, and, 
more important for our subsequent analysis, the value 770 = must be favored by the Etherington result. In other words, 
for a given data set, the likelihood of 770 must peak at 770 = to satisfy the cosmic relation. As we shall see, for those 
accepting the strict validity of the standard DD relation, our analysis suggests that galaxy clusters have an elliptical 
geometry. In principle, this kind of result is an interesting example of how a cosmological (global) condition correlates 
with the local physics. 

2. Galaxy cluster samples 

To constrain the values of 770, we consider two galaxy cluster samples for which DA(z) data were derived by combining 
their SZE and X-ray surface brightness observations. 
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Fig. 1. Galaxy clusters data. The open (blue) and filled (red) circles with the associated error bars represent, respectively, 
the Dc Filippis et al. (2005) and Bonamente et al. (2006) samples. 



The first sample is formed by 25 galaxy clusters compiled by De Filippis et al. (2005). By using an isothermal elliptical 
(3 model to describe the clusters, the DA(z) data was derived for two subsamples discussed in the literature. The first one, 
compiled by Reese et al. (2002), is a selection of 18 galaxy clusters distributed over the redshift interval 0.14 < z < 0.8. 
The second subsample of Mason et al. (2001) has 7 clusters from the X-ray limited flux sample of Ebeling et al. (1996). 
These 25 pieces of data are referred to as the elliptical sample. 

The second sample is defined by the 38 galaxy clusters observed by Bonamente et al. (2006), where the cluster 
plasma and dark matter distributions were analyzed assuming hydrostatic equilibrium model and spherical symmetry. 
This sample consists of clusters that have both X-ray data from the Chandra Observatory and SZE data from the 
BIMA/OVRO SZE imaging project, which uses the BerkeleyTllinois-Maryland Association (BIMA) and Owens Valley 
radio observatory (OVRO) interferometers to image the SZE. This dataset is termed the spherical sample. 

In Fig. [TJ we plot the elliptical and spherical galaxy cluster samples. Some authors have adopted these samples 
to estimate the galaxy cluster distances and measure the Hubble parameter by means of the SZE/X-ray technique 
(Bonamente et al. 2006; Cunha, Marassi & Lima 2007). However, since these samples are endowed with different geometric 
assumptions, our main interest here is to confront these underlying hypotheses with the validity of the DD relation. 

At present, there is no convincing evidence for deviations from the minimal cosmic concordance fiat model (Komatsu 
et al. 2011, Percival et al. 2010, Riess et al. 2009). In this case, the angular diameter distance is given by (Lima et al. 
2003; Cunha, Marassi & Lima 2007) 

„ , , n , 3000/1- 1 f z dz' 

where h — i?o/100 km s _1 Mpc -1 and the dimensionless function TL(z'; f2 m ) is given by 

H= [O m (l + z') 3 + (l-O m )] 1/2 . (5) 

In our statistical analysis (see next section), the Da(z) value was calculated by using the cosmological parameters given 
in the joint analysis carried out by Komatsu et al. (2011), namely ft\ = (1 — il m ) = 0.728 ± 0.015 and h = 0.704 ± 0.013. 



3. Analysis and Results 

We estimate the 770 parameter for each sample and both parameterizations of r)(z), namely, rj(z) = 1 + t/qz and T](z) = 
1 + r]oz/(l + z). To begin with, we evaluate the likelihood distribution function e~ x / 2 , where 

x 2 = y> W(z) - Vobsi*)] 2 (6) 

and T]g bs (z) = D A iata [z) / D a{z) (see Eq. ([3])). The statistical and systematic errors have been discussed by many authors 
(Mason et al, 2001; Reese et al. 2002, 2004; Bonamente et al. 2006). Statistical error contributions for galaxy clusters 
are: SZE point sources ±8%, X-ray background ±2%, Galactic N# < ±1%, ±15% for cluster asphericity, ±8% kinetic SZ 
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Fig. 2. The likelihood distribution functions for spherical (dotted-dashed and dashed red lines) and elliptical (dotted 
and solid blue lines) cases. See the text for values of 770. 



and for CMBR anisotropy < ±2%. On the other hand, the estimates of systematic effects are: SZ calibration ±8%, X-ray 
flux calibration ±5%, radio halos +3%, and X-ray temperature calibration ±7.5%. One may show that typical statistical 
errors can reach nearly 20%, in agreement with other works (Mason et al. 2001; Reese et al. 2002, 2004), whereas for 
systematics we also find typical errors around +12.4% and —12% (see also table 3 in Bonamente et al. 2006). 
In Fig. [51 we plot the likelihood distribution function for each galaxy cluster data. The basic results are: 

- For the elliptical sample, we obtain 770 = -0.056^;^ (xl.o.f. = °- 98 ) and Vo = -0.0881^4 (xl.o.f. = °- 97 ) in 68% C - L 
for a linear (blue solid line) and non linear (blue dotted line) parametrization, respectively. 

- For the spherical sample, we obtain 770 = -0.12^°;^ (xl.o.f. = °- 85 ) and Vo = -0.175±£o83 (Xd. .f. = °- 84 ) in 68% 
c.l. for a linear (red dashed line) and non linear (red dotted dashed line) parametrization. 

The uncertainties in rjo include (in quadrature) the systematic plus statistical errors. We can see that the D A lata (z) 
inferred based on the elliptical description agrees with the DD relation, whereas the other case, in which a spherical (3 
model was assumed to describe the clusters, is only marginally compatible with the DD relation. This result remains valid 
even when only clusters with z > 0.1 are considered in the elliptical sample. In this situation, we obtain ry = -0.044±°;} 

(Xd o f = 0.94) for the linear parametrization, and 770 = — 0.07^° (xd o f = 0.93) within la in the non-linear case. 

It thus follows that we have found no evidence of violation of the DD relation when the elliptical case is considered. 
However, the same kind of analysis shown in Fig. 2 points to a contradiction with the spherical symmetry hypothesis 
assumed in the Bonamente et al. (2006) sample. 



4. Conclusions 

We have explored some consequences of a deformed distance duality relation, r/(z) = Dl(1 + z)~ 2 /Da, based on obser- 
vations of Sunyaev-Zeldovich effect and X-ray from galaxy clusters. The consistency between the strict validity of the 
standard relation (r](z) = 1) and the assumptions regarding the geometry used to describe the galaxy clusters (elliptical 
and spherical f3 models) has been discussed. The rj(z) function was parametrized in two distinct forms, r\ — 1 + t/qz and 
77 = 1 + t]qz/(1 + z), where 770 is a constant parameter quantifying a possible departure from the strict validity of the 
duality relation. The basic idea pursued in this work is a simple one. The likelihood of the free parameter appearing in 
the proposed expressions for rj(z) should peak around 770 = when the distance duality relation is strictly obeyed. 

By comparing the De Filippis et al. (2005) (elliptical isothermal f3 model) and Bonamente et al. (2006) (spherical non- 
isothermal (3 model) samples with Da(z) obtained from ACDM (WMAP7), we show that the elliptical geometry is more 
consistent with no violation of the duality relation. The uncertainties in rjo included the systematic plus statistical errors 
from cluster data. In the case of an elliptical sample (see Fig. [2]), we found that 770 = — 0. 056j^i and 770 = —0.088^0^4 
for the linear and non-linear parametrization, respectively. However, the spherical sample (see Fig. 2) is only marginally 
compatible with rjo = — 0.12^Q Qgg and r/o = — 0.175^o'Qgg for linear and non-linear parametrization, respectively. Our 
analysis reveals that the elliptical model is compatible with the duality relation validity at la, whereas the spherical 
model is only marginally compatible at 3a. 

At this point, it is interesting to compare our results with those obtained by following a complementary approach 
(Holanda, Lima & Ribeiro 2010). The rj(z) function there was also parametrized as in the present work. However, the 
overall discussion was based on a model-independent cosmological test by considering Da{z) from galaxy clusters and 
the luminosity distances given by two sub-samples of SNe la taken from the constitution data (Hicken et al. 2009). Both 
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analyse are consistent with each other and suggest that the elliptical model is more compatible with the validity of the 
standard duality relation than the spherical case. 

Summarizing, the statistical analysis presented here provides additional evidence that the true geometry of clusters 
has an elliptical form. In principle, it is remarkable that a local property (the geometry of galaxy clusters) might be 
constrained by a global argument such as the one provided by the cosmological distance duality relation. In the near 
future, as more and larger data sets with smaller statistical and systematic uncertainties become available, the method 
proposed here (based on the validity of the distance duality relation) can improve the limits on the measurements of 
cluster geometries. 
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